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(57) A method of analysing electromagnetic survey data 
from an area of seafloor 6 that is thought or known 
to contain a subterranean hydrocarbon reservoir 12 
is described. The method includes providing 
electric field data and magnetic field data, for 
example magnetic flux density, obtained by at least 
one receiver 25 from a horizontal electric dipole 
(HED) transmitter 22 and determining a vertical 
gradient in the electric field data. The vertical 
gradient in the electric field data and the magnetic 
field data are then combined to generate combined 
response data. The combined response data is 
compared with background data specific to the area 
being surveyed to obtain difference data sensitive 
to the presence of a subterranean hydrocarbon 
reservoir. Because the combined response data are 
relatively insensitive to the transverse electric (TE) 
mode component of the transmitted signal, the 
method allows hydrocarbon reservoirs to be 
detected in shallow water where the TE mode 
component interacting with the air would otherwise 
dominate. Furthermore, because there is no mixing 
between the TE and transverse magnetic (TM) 
modes in the combined response data, data from 
all possible transmitter and receiver orientations 
may be used. 
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TITLE OF THE INVENTION 

ELECTROMAGNETIC S URVEYING FOR HYDROCARBON RESERVOrRS 
BACKGROUND OF THE INVENTION 

The invention relates to seafloor electromagnetic surveying for oil and other 
hydrocarbon reserves. 

Figure 1 schematically shows a surface vessel 14 undertaking controlled 
source electromagnetic (CSEM) surveying of a subterranean strata configuration using 
standard techniques [1]. The subterranean strata configuration in this example 
includes an overburden layer 8, an underburden layer 9 and a hydrocarbon reservoir 
12. The surface vessel 14 floats on the surface 2 of a body of water, in this case 
seawater 4 of depth h metres. A submersible vehicle 19 carrying a sourx;e in the form 
of a horizontal electric dipole HED transmitter 22 is attached to the surface vessel 14 
by an umbilical cable 16. This provides an electrical and mechanical connection 
between the submersible vehicle 19 and the surface vessel 14. The HED transmitter is 
supplied with a drive current so that it broadcasts an HED electromagnetic (EM) 
signal into the seawater 4. The HED transmitter is positioned a height z ' (typically 
around 50 metres) above the seafloor 6. The EM signals comprise transverse electric 
(TE) and transverse magnetic (TM) mode components. 

One or more remote receivers 25 are located on the seafloor 6. Each of the 
receivers 25 include an instrument package 26, a detector 24, a floatation device 28 
and a ballast weight (not shown). The detector 24 comprises an orthogonal pair of 
horizontal electric dipole detectors and an orthogonal pair of horizontal magnetic field 
detectors positioned a height z above the seafloor 6. The horizontal electric dipole 
detectors are sensitive to horizontal components of the electric fields induced by the 
HED transmitter in the vicinity of the receiver 25, and produce electric field detector 
signals therefrom. The horizontal magnetic field detectors are sensitive to horizontal 
components of the magnetic fields, for example the magnetic flux density, induced by 
the HED transmitter in the vicinity of the receiver 25, and produce magnetic field 



detector signals therefrom. The instrument package 26 records the detector signals for 
later analysis. Examples of suitable receivers are described by Constable [8] and US 5 
770 945 [9]. 

The HED transmitter 22 broadcasts EM signals that propagate outwards both 
into the overlying water column 4 and downwards into the seafloor 6 and the 
underlying strata 8, 9, 12. At practical frequencies for this method and given the 
typical resistivity of the respective media 4, 8, 9, 12, propagation occurs by diffusion 
of electromagnetic fields. The rate of decay in amplitude and the phase shift of the 
signal are controlled both by geometric spreading and by skin depth effects. Because 
in general the underlying strata 8, 9, 12 are more resistive than the seawater 4, skin 
depths in the underlying strata 8, 9, 12 are longer. As a result, electromagnetic fields 
measured by a receiver located at a suitable horizontal separation are dominated by 
those components of the transmitted EM signal which have propagated downwards 
through the seafloor 6, along within the underlying strata 8, 9, 12, and back up to the 
detector 24 rather than directly through the seawater 4. 

A sub-surface structure which includes a hydrocarbon reservoir, such as the 
one shown in Figure 1 , gives rise to a measurable increase in the horizontal electric 
field component amplitudes measured at the receiver relative to a sub-surface 
structure having only water-bearing sediments. This is because hydrocarbon reservoirs 
have relatively high resistivities (typically 100 fim) compared to other subterranean 
strata (typically 1 Qm) and so the EM signals are less attenuated. It is this 
enhancement in horizontal electric field amphtudes which has been used as a basis for 
detecting hydrocarbon reservoirs [1], 

It is important when surveying for hydrocarbon reservoirs to carefully consider 
the orientation of the current flows induced by a transmitted EM signal. The response 
of seawater and subterranean strata (which will typically comprise planar horizontal 
layers) to EM signals is generally very different for TE mode components of the 
transmitted signal, which excite predominantly horizontal current flows, and TM 
mode components, which excite significant components of vertical current flow. 



For TE mode components, the coupling between the layers comprising the 
subterranean strata is largely inductive. This means the presence of thin resistive 
layers (which are indicative of hydrocarbon reservoirs) does not significantly affect 
the EM fields detected at the surface as the large scale current flow pattern is not 
affected by the thin layer. On the other hand, for TM mode components, the coupling 
between layers includes a significant galvanic component (i.e. due to the direct 
transfer of charge between layers). For the TM mode even a thin resistive layer 
strongly affects the EM fields detected at the receiver since the large scale current 
flow pattern is interrupted by the resistive layer. It is known therefore that a significant 
component of the TM mode is required to satisfactorily perform an EM survey in the 
field of oil exploration. 

However, sole reliance on the sensitivity of the TM mode components to the 
presence of a thin resistive layer can lead to ambiguities. The effects on detected EM 
fields arising firom the presence a thin resistive layer can be indistinguishable fi-om the 
effects arising firom other realistic large scale subterranean strata configurations. In 
order to resolve these ambiguities it is known to determine the response of the 
subterranean strata to both TM mode components (i.e. inductively coupled) and TE 
mode components (i.e. galvanically coupled) [1]. 

The HED transmitter 22 shown in Figure 1 simultaneously generates both TE 
and TM mode components with the relative contribution of each mode to the signal at 
the receiver depending on the HED transmitter-receiver orientation. At receiver 
locations which are broadside to the HED transmitter axis, the TE mode dominates 
the response. At receiver locations which are inline with the HED transmitter axis, the 
TM mode is stronger (although the TE mode is still present) [1, 2, 3, 4]. The response 
at receiver locations in both the inline and broadside configurations is govemed by a 
combination of the TE and TM mode components, and these tend to work in 
opposition. 

Previous surveys [5, 6] have relied on this geometric splitting of the TE and 
TM mode components to determine the different response of the subterranean strata to 



the different modes. This is achieved by collecting electric field amplitude data for 
different transmitter-receiver alignments. This approach provides complementary 
horizontal electric field amplitude data sets which are differently sensitive to the TE 
and TM mode components of the transmitted EM signals. During analysis, these 
complementary data sets are combined to reveal differences between the TE mode and 
TM mode coupling between the transmitter and the receiver. These differences are 
indicative of the presence or not of a subterranean hydrocarbon reservoir. Because of 
the need to survey with multiple transmitter-receiver alignments, this approach 
requires a relatively large numbers of tow lines and receivers to ensure adequate 
coverage. 

Figure 2 shows in plan view an example survey geometry for collecting 
horizontal electric field component data to be analysed according to known methods. 
Sixteen receivers 25 are laid out in a square grid on a section of seafloor 6 above a 
subterranean reservoir 56 having a boundary indicated by a heavy line 58. The 
orientation of the subterranean reservoir is indicated by the cardinal compass points 
(marked N, E, S and W for North, East, South and West respectively) marked in the 
upper right of the figure. To perform a survey, a transmitter starts firom location 'A' 
and is towed along a path indicated by the broken hne 60 through location 'B' until it 
reaches location 'C which marks the end of the survey path. As is evident, the tow 
path first covers four parallel paths aUgned with the North-South direction to drive 
over the four "columns" of the receivers. This part of the survey path moves from 
location 'A' to 'B'. Starting firom location 'B', the survey path then covers four paths 
aligned with the East- West direction which drive over the four "rows" of receivers. 
Each receiver is thus driven over in two orthogonal directions. The survey is 
completed when the transmitter reaches the location marked *C\ 

During the towing process, each of the receivers 25 presents several different 
orientation geometries with respect to the transmitter. For example, when the 
transmitter is directly above the receiver position Dl and on the North-South aligned 
section of the tow path, the receivers at positions D5, D6 and D7 are at different 



separations in an inline position (i.e. aligned with the dipole axis of the HED 
transmitter), the receivers at positions D2, D3 and D4 are at different horizontal 
separations in a broadside position and the receiver at positions D8 and D9 are in- 
between. However, when the transmitter later passes over the receiver position Dl 
when on the East-West aligned section of the tow path, the receivers at positions D5, 
D6 and D7 are now in a broadside position, and the receivers at position D2, D3 and 
D4 are in an inline position. Thus, in the course of a survey, and in conjunction with 
the positional information of the transmitter, data from the receivers can be used to 
provide details of the signal transmission through the subterranean strata for a 
comprehensive range of distances and orientations between transmitter and receiver, 
each with varying TM mode and TE mode contributions to the signal propagation. 

In addition to requiring relatively complex tow paths, another problem with 
known survey and analysis techniques is they do not provide good results for surveys 
made in shallow waters. This is due to the presence of an *airwave' component in the 
EM fields induced by the HED transmitter at the receiver. This airwave component is 
due to EM signals from the HED transmitter which interact with the air. Since air is 
non-conducting and hence causes little attenuation, the airwave component can 
dominate the fields at the receiver. The airwave component is principally due to the 
TE mode components. This is because the TE mode components are efficiently 
inductively coupled across the seawater-to-air interface. The TM mode components, 
on the other hand, do not couple well across this boimdary and consequently do not 
contribute significantly to the airwave component. The airwave component contains 
little information about subterranean resistivity. Accordingly, if the airwave 
contributes a significant component to the EM fields induced by the HED transmitter 
at the receiver, the sensitivity of the technique to subterranean resistivity structures, 
such as hydrocarbon reservoirs, is greatly reduced. The path of an example airwave 
component is shown in Figure 1 by a dotted line labelled AW. The magnitude of the 
airwave component is reduced only by geometric spreading. This is because air is 
non-conducting. However, as with other components, the airwave component is 
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strongly attenuated by its passage through the seawater. This means that in relatively 
deep water (large h) the airwave component is not very significant at the receiver and 
as such does not present a major problem. However in shallow water (small h) the 
airwave component does not pass through as much seawater and thus makes a larger 
5 contribution to the EM fields induced by the HED transmitter at the receiver. This 
contribution becomes greater still at increasing transmitter-receiver horizontal 
separations. This is because (other than due to geometric spreading) the strength of the 
airwave component is relatively constant over a wide range of horizontal separations 
since any extra distance travelled by the airwave component is almost exclusively in 
1 0 the non-attenuating air. Other components of the EM fields induced by the HED at the 
receiver, such as those which pass through the subterranean strata and are of interest, 
travel through lower resistivity media and become increasing attenuated as they travel 
further. For these reasons, the airwave component tends to dominate the EM fields 
induced by the HED transmitter at the receiver for surveys made in shallow water, 
1 5 especially at long transmitter-receiver horizontal separations. 

The existence of the airwave as a dominant component of the detector signals 
limits the applicability of the above described surveying and analysis techniques. In 
shallow water the transmitter-receiver separations over which the techniques can be 
applied is much reduced. This not only leads to a need to employ more receiver 
20 locations to adequately cover a given area, but also limits the depth beneath the 
seafloor to which the technique is sensitive. This can mean that a buried hydrocarbon 
reservoir in shallow water may not be detectable, even though the same reservoir 
would be detected in deeper water. 

Figure 3A is a graph schematically showing results of one-dimensional 
25 modelling of two example EM surveys of the kind shown in Figure 1. One example 
corresponds to a survey performed in deep water (dotted line) and the other to a 
survey performed in shallow water (solid line). For each model survey the amplitude 
of an electric field component induced at the receiver in response to the HED EM 
transmitter is calculated per unit transmitter dipole moment and is plotted as a 



function of horizontal separation r between the HED transmitter and the receiver. For 
both model surveys, the subterranean strata configuration is a semi-infinite 
homogeneous half space of resistivity 1 Qm, In the deep-water example, the 
subterranean strata configuration is located beneath an infinite extent of seawater. hi 
the shallow-water example, it is located beneath a 500-meter depth of seawater. hi 
both cases the seawater has resistivity 0.3 Qm. The transmitter and receiver are 
separated along a Une which runs through the axis of the HED transmitter (inline 
orientation). It is the component of detected electric field resolved along this direction 
which is plotted in Figure 3 A. The HED transmitter is driven by an altemating current 
(AC) drive signal at a frequency of 0.25 Hz. 

The effect of the airwave component on the amplitude of EM fields induced by 
the HED transmitter at the receiver is clear, hi the deep-water model survey, where 
there is no airwave component, the calculated electric field amplitude falls steadily 
with increasing horizontal separation. Li the shallow-water model, however, where 
there is a strong airwave component, the rate of amplitude reduction sharply decreases 
at a transmitter-receiver horizontal separation of about 5000 m. Figure 3B is a plot 
showing the ratio, p, of the two curves shown in Figure 3A. The large deviations from 
unity seen in Figure 3B highlight the difference between these curves. Since the only 
difference between the two model surveys is the presence or not of an airwave 
component, the ratio plotted in Figure 3 A effectively shows the relative strength of the 
airwave component in the detected signal compared to that which passes through the 
subterranean strata for the shallow- water model survey. 

It is apparent from Figures 3A and 3B that at all but the very shortest 
horizontal separations the detected electric field is significantly larger in the shallow- 
water model. For example, at a horizontal separation of 2500 m, the amplitude of the 
detected signal in the deep-water model survey is around 10'^^ V/Am^. In the shallow- 
water model survey it is higher at around 10"*' ^ V/Am^. This is due to the additional 
contribution of the airwave component. This level of increase shows that the airwave 
component has an amplitude more than double that of the component which has 



passed through the subterranean strata, and accordingly over two-thirds of the detector 
signal carries almost no information about the subterranean strata. At greater 
horizontal separations the airwave component dominates even more. In particular, it 
becomes especially pronounced beyond around 5000 m. At this point there is a break 
in the rate at which the detected electric field amplitude falls with increasing 
horizontal separation. At a horizontal separation of around 7000 m, the airwave 
component in the shallow-water example has an amplitude around twenty times 
greater than that of the signal which passes through the subterranean strata. This 
clearly imposes high requirements for the signal-to-noise ratio of data collected over 
these sorts of horizontal separations, as is generally the case when a small signal rides 
on a large background. It is apparent that the airwave significantly limits the 
usefulness of these surveying and analysis techniques in shallow water. 

While this survey method has been demonstrated to provide good results in 
practice, as noted above some limitations have been identified. 

Firstly, since the TE and TM mode components cannot be easily separated 
there will generally be a level of cross-talk between them at a receiver. This may lead 
to ambiguities in the results. 

Secondly, in order to obtain survey data from both inline and broadside 
geometries, the HED transmitter needs to be re-oriented at each HED transmitter 
survey location. This requires the surface vessel to make multiple passes over 
broadcast locations and can lead to complex and long tow patterns. 



SUMMARY OF THE INVENTION 



According to a first aspect of the invention, there is provided a method of 
analysing results from an electromagnetic survey of an area that is thought or known 
to contain a subterranean hydrocarbon reservoir, comprising: providing electric field 
data and magnetic field data obtained by at least one receiver fi-om at least one 
horizontal electric dipole (HED) transmitter; determining a vertical gradient in the 
electric field data; and combining the vertical gradient in the electric field data with 
the magnetic field data to generate combined response data. 

In this case, references to vertical and horizontal indicate that a significant 
component, preferably a major component, of the respective signals should be aligned 
with the vertical and horizontal axes. It is not necessary that the signals are perfectly 
aligned to the vertical and horizontal axes, although close alignment is preferred to 
provide a strong signal and reduce the complexity of analysis, for example alignment 
within H-/-30^ is desirable. 

By generating combined response data in the analysis, it is possible to analyse 
survey results taken in shallower water than has previously been possible. This is 
because the combined response data are not sensitive to transverse electric (TE) mode 
components which have propagated through air and which tend to dominate survey 
results analysed with previous methods. In addition, because suitable combined 
response data can be formed for all transmitter-receiver orientations, a complete 
survey can be carried out witii a single tow. In contrast, with previously proposed EM 
survey analysis methods based on use of HED transmitters, two sets of survey data 
had to be obtained with different tows over different paths in order to probe the 
response of an area of interest galvanically and inductively. 

Combined response data of this kind are fimctionally similar to vertical 
electric field data. Accordingly, the combined response data provide similar benefits 
to those achievable with data fi-om vertical electric field detectors, such as the benefits 
described in the inventors' co-pending UK patent application number 0313376.6 [10]. 
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With the present invention, these are achieved without relying on vertical field 
measurements. This is beneficial because vertical field measurements can be more 
susceptible to noise, particularly motion-induced noise caused by undersea water 
currents. 

The electric field data may include a horizontal component of electric field 
resolved along a first direction and the magnetic field data may include a horizontal 
component of magnetic flux density, resolved along a second, different, direction. The 
first and second directions may be orthogonal. Furthermore, in some examples the 
first direction may be parallel to a line connecting the HED transmitter to the receiver 
while in other examples it may be perpendicular to such a line. 

Although specialist equipment can be used to collect the data, suitable data can 
also be collected using existing equipment. This means the method may be applied to 
existing data. The analysis method allows for relatively simple combinations, for 
example, linear combinations, of the vertical gradient in the electric field data and the 
magnetic field data to be made to generate appropriate combined response data. 

Vertical gradients in the electric field may be determined in several ways. One 
way is by comparing electric field data for different heights. For example for different 
heights of transmitter and/or detector. Two or more different heights may be used. 
Another way is to compare the electric field data with data simulated using a 
backgroimd model. The data simulated using a background model may include a 
boundary condition, such as an assumed or modelled electric field at a boundary 
between the seawater and air overlying the area surveyed, for example. Yet another 
way is to determine the vertical electric field gradient at or in the vicinity of a first 
receiver by comparing electric field data from the first receiver when the transmitter is 
above a second receiver with electric field data fl-om the second receiver when the 
transmitter is above the first receiver, and applying a predetermined adjustment to the 
electric field data fi'om second receiver. 

The method may further comprise providing background data specific to the 
area being surveyed; and comparing the combined response data with the background 
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data to obtain difference data sensitive to the presence of a subterranean hydrocarbon 
reservoir. 

This can be beneficial since the comparison of the combined response data 
with background data can help to determine whether features of the combined 
response data are indicative of a hydrocarbon reservoir or arise as a result of some 
other local background structure configuration. Background data may be obtained by 
modelling the EM survey performed to obtain the combined response data with a 
model backgrovind subterranean strata configuration. The background model strata 
configuration should preferably be a close match to the actual background stmcture in 
the area being surveyed. 

The background data may be obtained in several ways, for example from a 
controlled source electromagnetic survey, from a magneto-telluric electromagnetic 
siirvey, from another similar survey taken at a different time, or from a rock formation 
model. If a rock formation model is used it should preferably include resistivity, and 
may be derived from a combination of geological data and resistivity data. The 
geological data can be from seismological surveying and the resistivity data from well 
logging. Other sources of inforaiation, such as neutron data or other porosity estimates 
from well logs, could also be used. 

In some examples, the background data may be obtained from the electric and 
magnetic field data used to provide the combined response data. This can be achieved 
by combining the data in a different way. For example, determining a vertical gradient 
in the magnetic field data and combining it with electric field data. The vertical 
gradient in the magnetic field data can be determined in ways similar to those 
described above for the vertical gradient in the electric field data. 

The difference data may represent the difference between the combined 
response data and the background data as a function of position within the area 
surveyed, and the analysis may include identifying a location of a boundary of a 
subterranean hydrocarbon reservoir. 
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According to a second aspect of the invention there is provided a computer 
program product bearing machine readable instructions for implementing a method of 
analysing results from an electromagnetic survey according to the first aspect of the 
invention. 

5 According to a third aspect of the invention there is provided a computer 

apparatus loaded with machine readable instructions for implementing the method of 
analysing results from an electromagnetic survey according to the first aspect of the 
invention. 

According to a fourth aspect of the invention there is provided a method of 
10 planning an electromagnetic survey of an area that is thought or known to contain a 
sr|)terranean hydrocarbon reservoir, comprising: creating a model of the area to be 
surveyed including a rock formation containing a postulated hydrocarbon reservoir, 
and a body of water above the rock formation; setting values for water depth, depth of 
the postulated hydrocarbon reservoir, and resistivity structure of the rock formation; 
15 and performing a simulation of an electromagnetic survey in the model of the survey 
area by calculating electric field data and magnetic field data obtained by at least one 
simulated receiver detecting signals from a simulated horizontal electric dipole (HED) 
transmitter; determining a vertical gradient in the electric field data; and combining 
the vertical gradient in the electric field data with the magnetic field data to generate 
20 combined response data. 

The method may further comprise adjusting the model to remove the 
postulated hydrocarbon reservoir and repeating the simulation to obtain background 
data for comparison with the combined response data. 

Repeated simulations for a number of transmitter-receiver horizontal 
25 separations and frequencies of signal can be performed in order to allow optimum 
surveying conditions in terms of transmitter-to-receiver distance and frequency of EM 
signal for probing the hydrocarbon reservoir to be selected when performing an 
electromagnetic survey. The effects and usefulness of differing receiver array 
configurations and transmitter tow paths can also be modelled. 
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According to a fifth aspect of the invention there is provided a computer 
program product bearing machine readable instructions for implementing the method 
of planning an electromagnetic survey according to the fourth aspect of the invention. 

According to a sixth aspect of the invention there is provided a computer 
5 apparatus loaded with machine readable instructions for implementing the method of 
planning an electromagnetic survey according to the fourth aspect of the invention. 

According to a seventh aspect of the invention there is provided an 
electromagnetic survey method applied to a survey area that is thought or known to 
contain a subterranean hydrocarbon reservoir, comprising: providing at least one 
10 transmitter and at least one detector for transmission and detection of electromagnetic 
signals; and obtaining data with transmission and/or detection at a plurality of 
different heights over the survey area, so that the data allow comparison of 
electromagnetic signals transmitted and^or received at different vertical displacements. 

Such a survey method provides data which allow gradients in electric field 
1 5 data to be determined such that the data may be analysed according to the methods of 
the first aspect of the invention. 

The transmission and/or detection at a plurality of different heights may 
comprise detection at a plurality of different heights. The detection at a plurality of 
different heights may be made simultaneously using multiple detectors. 
20 Because of the reciprocity between transmitter and detector, the transmission 

and/or detection at a plurality of different heights may equally comprise transmission 
at a plurality of different heights. The transmission at a plurality of different heights 
may made simultaneously by multiple transmitters, or at different times by a single 
transmitter. For example, a mobile transmitter may be first towed along a tow path at 
25 one height and then re-towed along at least parts of the tow path at a different height. 

According to an eighth aspect of the invention there is provided an 
electromagnetic receiver for use in an electromagnetic survey of an area that is 
thought or known to contain a subterranean hydrocarbon reservoir, the receiver, when 
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normally deployed, being operable to measure electric fields at two or more different 
heights such that a vertical gradient in electric field may subsequently be determined. 

Such a receiver is suitable for performing a survey according to the seventh 
aspect of the invention. 

The receiver may be further operable to measure magnetic fields at two or 
more different heights such that a vertical gradient in magnetic field may subsequently 
be determined. 

The receiver, when normally deployed, may comprise first and second pairs of 
horizontally aligned electric dipole detectors extending in different horizontal 
directions, the pairs being positioned at different heights. 

The receiver, when normally deployed, may fiuther comprise first and second 
pairs of magnetic field detectors for detecting magnetic fields in different horizontal 
directions, the pairs being positioned at different heights. 

According to a ninth aspect of the invention there is provided and 
electromagnetic source for use in an electromagnetic survey of an area that is thought 
or known to contain a subterranean hydrocarbon reservoir, the source comprising two 
transmitters which in normal use are at different heights. 

The source may be a mobile source, for example a pair of HED transmitters 
towed by a single submersible vehicle, or a pair of submersible vehicles each toeing a 
single transmitter. 

Source according to the ninth aspect of the invention are suitable for 
performing a survey according to the seventh aspect of the invention. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

For a better understanding of the invention and to show how the same may be 
carried into effect reference is now made by way of example to the accompanying 
drawings, in which: 

Figure 1 shows in schematic vertical section a surface vessel undertaking an 
EM survey in deep water according to standard techniques; 

Figure 2 is a schematic plan view showing an example survey geometry for 
providing data to be analysed according to a previously proposed method in which 
sixteen receivers are laid out on a section of seafloor above a subterranean reservoir; 

Figure 3A is a graph plotting detector signals calculated from two model 
surveys analysed according to the previously proposed method, one performed in deep 
water (dotted line) and one performed in shallow water (soUd line); 

Figure 3B is a graph plotting the ratio of the two curves shown in Figure 3 A; 

Figure 4 is a plan view detailing a cylindrical polar co-ordinate system; 

Figures 5A-5F show equations (Equations 1 to 6) corresponding to solutions to 
Maxwell's equations for an HED transmitter in a finite seawater layer overlying a one- 
dimensional subterranean resistivity structure for the radial (r), azimuthal ((])) and 
vertical (z) components of the electric (E) and magnetic (B) fields; 

Figure 5G shows an equation (Equation 7) defining a linear combination of a 
vertical gradient in radial electric field data (Er) and azimuthal magnetic field data 
(B^) used in a method of analysing survey data according to a first embodiment of the 
invention; 

Figure 5H shows an equation (Equation 8) defining a linear combination of a 
vertical gradient in the azimuthal electric field data (E^j^) and the radial magnetic field 
data (Br) used in a method of analysing survey data according to a second embodiment 
of the invention; 

Figure 6A schematically shows in vertical section a surface vessel undertaking 
an EM survey according to an embodiment of the invention; 
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Figure 6B schematically shows in vertical section a surface vessel undertaking 
an EM survey according to another embodiment of the invention; 

Figure 7A schematically shows in vertical section a survey using similar 
equipment to that shown in Figure 1 ; 

Figure 7B schematically shows in vertical section the survey shown in Figure 
7A at a later time; 

Figure 8 shows in schematic vertical section a model uniform background 
subterranean strata configuration; 

Figures 9A and 9B show the modelled amplitude and phase of the radial 
component of electric field data collected during an EM survey of the model 
subterranean strata configuration shown in Figure 8 for a range of water depths; 

Figures lOA and lOB show the modelled ampUtude and phase of a 
combination of electric and magnetic field data collected during an EM survey which 
provides inline TM mode decomposition data according to a first embodiment of the 
invention for a range of water depths; 

Figures IOC and lOD show the modelled amplitude and phase of a 
combination of electric and magnetic field data collected during an EM survey which 
provides broadside TM mode decomposition data according to a second embodiment 
of the invention for a range of water depths; 

Figure 1 1 shows in schematic vertical section a model hydrocarbon-reservoir 
bearing subterranean strata configuration; 

Figures 12A and 12B show the modelled amphtude and phase of the radial 
component of electric field data collected during an EM survey of the model 
subterranean strata configuration shown in Figure 11 for a range of hydrocarbon- 
reservoir resistivities; 

Figures 13A and 13B show the modelled amplitude and phase of the inline 
TM mode decomposition data for an EM survey of the model subterranean strata 
configuration shown in Figure 11 for a range of hydrocarbon-reservoir resistivities; 
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Figures 13C and 13D show the modelled amplitude and phase of the broadside 
TM mode decomposition data for an EM survey of the model subterranean strata 
configuration shown in Figiure 11 for a range of hydrocarbon-reservoir resistivities; 

Figure 14A is a graph schematically showing the ratio of the curves plotted in 
Figures 12 A (radial electric field component), 13A (inline TM mode decomposition) 
and 13C (broadside TM mode decomposition) for a hydrocarbon reservoir of 
resistivity p = 100 Clm to the corresponding curves for which there is no detectable 
hydrocarbon reservoir; 

Figure 14B is a graph schematically showing the difference in phase between 
the curves plotted in Figures 12B (radial electric field component), 13B (inline TM 
mode decomposition) and 13D (broadside TM mode decomposition) for a 
hydrocarbon reservoir of resistivity p = 1 00 Qm and the respectively corresponding 
curves for which there is no detectable hydrocarbon reservoir; 

Figure 14C is a graph schematically showing the same curves as shown in 
Figure 14A, but on a magnified scale; 

Figure 15A shows in schematic vertical section another model hydrocarbon- 
reservoir bearing subterranean strata configuration; 

Figure 15B shows in schematic vertical section a model increasing-resistivity 
subterranean strata configuration in which resistivity progressively increases with 
depth; 

Figure 1 6 is a graph schematically showing the ratios of the inline TM mode 
decomposition and the broadside TM mode decomposition calculated for both the 
model subterranean strata configurations shown in Figures 15A and 15B to those 
calculated for the background model subterranean strata configuration shown in 
Figure 8 for infinite water depth; 

Figure 17A shows an equation (Equation 9) defining a linear combination of 
radial electric field data (Er) and a vertical gradient in azimuthal magnetic field data 
(B4.); 
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Figure 1 7B shows an equation (Equation 1 0) defining a linear combination of 
azimuthal electric field data (E^) and a vertical gradient in radial magnetic field data 

(Br); 

Figure 1 8 A is a graph schematically showing the ratios of an inline TE mode 
decomposition and a broadside TE mode decomposition calculated for both the model 
subterranean strata configurations shown in Figures 15A and 15B to those calculated 
for the background model subterranean strata configuration shown in Figure 8 for 
infinite water depth; 

Figure 18B is a graph schematically showing the same curves as shown in 
Figure 1 8 A, but on a magnified scale; and 

Figure 1 9 is a schematic plan view showing an example survey geometry for 
obtaining electric and magnetic field data to be analysed according to embodiments of 
the invention in which sixteen receivers are laid out on a section of seafloor above a 
subterranean reservoir. 
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DETAILED DESCRIPTION 

Figure 4 is a schematic plan view showing a co-ordinate system for describing 
the relative placement of an HED transmitter 22 and a receiver 25 of the kind shown 
5 in Figure 1 . The position of the receiver 25 with respect to the HED transmitter 22 is 
most suitably described in cylindrical polar co-ordinates, with the centre of the HED 
transmitter 22 providing the origin of the co-ordinate system. The position of the 
receiver 25 is defined by an azimuthal angle <|> and a distance r. The angle ^ is 
measured clockwise from a line passing through, and running parallel to, the HED 

10 transmitter axis, as indicated in Figure 4 by the line marked <|) = 0°. A receiver placed 
along this line, i.e. such that is has an azimuthal angle ^ of 0^, is referred to as being 
in an inline position. A receiver with an azimuthal angle ^ of 90^, such that it lies on 
the line marked ^ = 90® in Figure 4, is referred to as being in a broadside position. The 
electric field strength at the receiver measured by the orthogonal pair of horizontal 

15 electric dipole detectors is resolved into a radial component Er and an orthogonal 
azimuthal component E<^, as indicated in the figure. The magnetic flux density at the 
receiver measured by the orthogonal pair of horizontal magnetic field detectors is 
similarly resolved into a radial component Br and an orthogonal azimuthal component 
B^. The axial co-ordinate z extends vertically away from the seafloor. 

20 The fundamental equations governing electromagnetic induction in the earth 

are Maxwell's equations. At frequencies typically used in CSEM surveys 
displacement currents can be neglected to give: V.5 = 0, V xE-^i(oB = 0 and 
Vx£-//q<t^ = //q/, where E is the electric field strength, B is the magnetic flux 

density, a is the conductivity of the medium, no is the magnetic permeabiUty which is 
25 assumed to take its free space value, J is the source current density, and a single 
Fourier component proportional to e'^ is considered. Maxwell's equations can be 
solved numerically in two- or three-dimensions for a point HED transmitter, however 
a closed form exists only for one-dimensional structures. Chave & Cox [7] derive a 
solution for the case of an HED transmitter in an infinite depth of seawater for a one- 
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dimensional subterranean strata configuration (i.e. in which resistivity varies only in 
the vertical z-direction). 

The inventors have performed an extension of the analysis presented in Chave 
& Cox [7] to model an HED transmitter in a finite depth h of seawater. Solving 
Maxwell's equations for an HED transmitter in a finite seawater layer overlying a one- 
dimensional subterranean resistivity structure provides equations for the radial (r), 
azimuthal ((j)) and vertical (z) components of the electric field (E) and magnetic flux 
density (B) as shown in Figures 5A to 5F. Although this modelling has been 
performed for a one-dimensional strata configuration, similar modelling may be 
performed in two- or three- dimensions. 

Where the " ± " or the " + " optional operator appears in the equations shown 
in Figures 5 A to 5F, the upper symbol is used when z' > z and the lower symbol 
when z ' < z. In these equations, z ' and z are the heights of the HED transmitter and 
detector above the seafloor respectively, h is the depth of the seawater, jiq is the 
permeability of firee space, P is the transmitter dipole moment, 

-toiler) = ^,,^=^(krr' and Mkr) = ^^^J0-^ikr)-- are zeroth and 
first order Bessel functions respectively, p© is the resistivity of the seawater, k is a 

parameter analogous to the wave number in a Fourier integral, fi^ = W -^^^ , 

V Po 

Rl™ and R are coefficients defining the TM and TE mode interaction with the 
seafloor which depend on the resistivity structure of the subterranean strata 
configuration, and Ra is a coefficient defining the TE mode interaction with the air. 

In the presentation of Equations 1, 2, 4 and 5 (which are the equations 
describing the horizontal components of the fields) in Figures 5A, 5B, 5D and 5E, the 
equations are shown split over four lines of text with each line of text having a left 
and a right component. The left component on each line is marked "TM" and results 
from the TM mode component of the transmitted signal and the right component is 
marked "TE" and results from the TE mode component of the transmitted signal. 
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As previously noted, the airwave component is principally due to the 
interaction of the TE mode with the air, i.e. determined by the Ra^^ coefficient. As can 
be seen form Equations 1 and 2, Er and E^^ include both TM and TE components and 
so are affected by the airwave. This is why known methods of analysing results firom 
5 CSEM surveys which are based on electric field amplitude enhancement do not work 
well in shallow water. 

Equation 7, shown in Figure 5G, defines a linear combination of a vertical 
gradient in the radial electric field data Er and the azimuthal magnetic field data 
used in a method of analysing results according to a first embodiment of the invention. 
10 Combinations of the electric and magnetic field data such as this are referred to as 
combined response data. 

Equation 8, shown in Figure 5H, defines another linear combination of the 
vertical gradient in the azimuthal electric field data E4, and the radial magnetic field 
data Br used in a method of analysing results according to a second embodiment of the 
1 5 invention. 

One way of obtaining the vertical gradients in the electric field data is by 
sampling the electric field at two (or more) different heights during surveying, for 
example by deploying multiple receivers at different heights or a single receiver 
having multiple detectors at different heights. 

20 Figure 6A schematically shows a surface vessel 14 undertaking controlled 

source electromagnetic (CSEM) surveying of a subterranean strata configuration using 
a survey method according to one embodiment of the invention. The surface vessel 14 
floats on the surface 2 of a body of water, in this case seawater 4 of depth h metres. A 
submersible vehicle 19 carrying a source in the form of an HED transmitter 22 is 

25 attached to the surface vessel 14 by an umbilical cable 16 providing an electrical and 
mechanical connection between the submersible vehicle 19 and the surface vessel 14. 
The HED transmitter is supplied with a drive current so that it broadcasts an HED EM 
signal into the seawater 4. The HED transmitter is positioned a height z ' (typically 
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around 50 metres) above the seafloor 6. The surface vessel 14, submarine 19, 
umbilical 16 and HED transmitter 22 may be conventional. 

One or more remote receivers 125 are located on the seafloor 6. Each of the 
receivers 25 includes an instrument package 126, a lower detector 124, an upper 
detector 129, a floatation device 128 and a ballast weight (not shown). The upper and 
lower detectors 124, 129 each comprises an orthogonal pair of horizontal electric 
dipole detectors and an orthogonal pair of horizontal magnetic field detectors. The 
horizontal electric dipole detectors are sensitive to horizontal components of the 
electric fields induced by the HED transmitter in the vicinity of the receiver 125, and 
produce electric field detector signals therefrom. The horizontal magnetic field 
detectors are sensitive to horizontal components of the magnetic fields, for example 
the magnetic flux density, induced by the HED transmitter in the vicinity of the 
receiver 125, and produce magnetic field detector signals therefirom. The instrument 
package 126 records the detector signals for later analysis. The lower detector is 
positioned a height z\ above the seafloor 6 and the upper detector is positioned at a 
height Zu. Typically the difference between zi ad Zu will be around 1 0 meters, though 
other separations may also be used. A comparison of signals measured by the upper 
and lower detectors 129, 124 allows a gradient in electric field data to be detennined 
based on their vertical separation. Although in this example, the upper and lower 
detectors 124, 129 are substantially identical. It will be appreciated that in other 
examples they may be different, for example, one or other of them might not include 
magnetic field detectors. In some examples detectors at more than two height may be 
used. 

Because of the reciprocity between receiver and transmitter, vertical gradients 
in electric field data can also be determined firom data collected by a receiver having a 
detector positioned at only a single height, but firom a source having transmitters 
positioned at different heights, for example two or more different heights. 

Figure 6B schematically shows a surface vessel 14 undertaking controlled 
source electromagnetic (CSEM) surveying of a subterranean strata configuration using 
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a survey method according to another embodiment of the invention. The surface 
vessel 14 floats on the surface 2 of a body of water, in this case seawater 4 of depth h 
metres. A submersible vehicle 19 tows an electromagnetic source comprising an upper 
HED transmitter 122 and a lower HED transmitter 123. The submarine is attached to 
the surface vessel 14 by an umbilical cable 16 providing an electrical and mechanical 
connection between the submersible vehicle 19 and the surface vessel 14. The HED 
transmitters comprising the source are supplied with drive currents so that they 
broadcasts HED electromagnetic signals into the seawater 4. The upper HED 
transmitter 122 is positioned a height Zu' (typically around 50 metres) above the 
seafloor 6. The lower HED transmitter 123 is positioned a height Z\ ' (typically around 
10 metres lower than the upper HED transmitter 122). In other example, the upper and 
lower transmitters may be towed by separate submersible vehicles. 

One or more remote receivers 25 are located on the seafloor 6. Each of the 
receivers 25 include an instrument package 26, a detector 24, a floatation device 28 
and a ballast weight (not shown). The detector 24 comprises an orthogonal pair of 
horizontal electric dipole detectors and an orthogonal pair of horizontal magnetic field 
detectors positioned a height z above the seafloor 6. The horizontal electric dipole 
detectors are sensitive to horizontal components of the electric fields induced by the 
HED transmitter in the vicinity of the receiver 25, and produce electric field detector 
signals therefirom. The horizontal magnetic field detectors are sensitive to horizontal 
components of the magnetic fields, for example the magnetic flux density, induced by 
the HED transmitter in the vicinity of the receiver 25, and produce magnetic field 
detector signals therefrom. The instrument package 26 records the detector signals for 
later analysis. 

During surveying, the upper and lower HED transmitters may be driven 
alternately at the same frequency, or simultaneously at different frequencies such that 
the instrument package can distinguish signals received from each one. A comparison 
of signals measured from the upper and lower HED transmitters 122, 123 allows a 
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gradient in electric field data to be determined based on the vertical separation of the 
transmitters. 

In other examples, a vertical gradient in the electric field data may be obtained 
by using conventional receivers and a conventional HED transmitter which is towed 
along a path at one height, for example 50 metres above the seafloor, and then 
subsequently towed so as to re-trace at least portions of path at a different height, for 
example around 1 0 metres higher or lower. 

It is also possible to use standard survey apparatus which provides electric and 
magnetic field data at a single height z at, or near, the seafloor for a single HED 
transmitter at a height z ' above the seafloor. One benefit of this approach is that data 
previously collected using conventional survey methods may be analysed using 
methods according to the present invention. One way to determine vertical gradients 
in the electric field using such data is to compare the data with data simulated using a 
background model. This can be done, for example, to provide a boimdary condition 
which may be applied to the data. For example, the electric fields at the seawater-air 
boundary interface may be modelled and Srom this, and the depth of the seawater, an 
estimate of the vertical gradient in the electric field may be determined. Altematively, 
gradients in the electric field may be inferred by comparing data from a pair of 
receivers as the HED transmitter moves from being positioned above one receiver to 
being positioned above the other. 

Figure 7 A schematically shows in vertical section a CSEM survey made using 
equipment similar to that shown in Figure 1. There are two receivers 25, one located 
at position B and one located at position D. The detectors for each receiver are at a 
height z above the seafloor 6 (i.e. z is the effective height at which the electric field 
and magnetic flux density are measured). The HED transmitter 22 is at a position A. 
Position A is directly above position B at a height z' above the seafloor 6. Also 
marked is a position C directly above position D, again at a height z' above the 
seafloor 6. The receivers at locations B and D are separated by a distance r. In this 
configuration, we consider the electric fields and magnetic flux densities measured by 
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the receiver at D (i.e. at a height z above the seafloor) in response to the broadcast 
signal from the HED transmitter at A. 

Figure 7B is similar to Figure 7A but corresponds to a later time, the HED 
transmitter having moved from location A to location C during its tow. In this 
configuration, we consider the electric fields and magnetic flux densities measured by 
the receiver at B (again at a height z above the seafloor) in response to the broadcast 
signal from the HED transmitter at C. Because of the reciprocity between transmitter 
and receiver, this configuration is equivalent to a configuration in which the 
transmitter and receiver positions are swapped. That is to say. Figure 7B corresponds 
to a transmitter at location B and a receiver at location C. Accordingly, this 
configuration provides data equivalent to having a receiver at a height z ' above the 
seafloor and a transmitter at a height z above the sea floor over the same horizontal 
range r as shown in Figure 7A. Using a far-field approximation (i.e. where skin depth 
dominate in the seawater), the transmitter at B (height z) can be effectively converted 
into a transmitter at A (height z') by scaling the measured field amplitudes and 
shifting their phase by appropriate amounts. For example, to a first approximation an 
amplitude scaling factor and shift in phase calculated using tiie EM signal skin deptii 
in the seawater or in a medium with an apparent resistivity calculated &om the 
measured of electric and magnetic field, may be used. Following this adjustinent, the 
data seen for ttie configuration shown in Figure 7B is equivalent to tiiat which would 
be seen with the HED ti-ansmitter at location A and tiie receiver at location C. This 
data can then be compared with data collected for tiie configuration shown in Figure 
7A (i.e. transmitter at location A, receiver at location D) to provide the necessary 
gradients in tiie fields over tiie height range z ' to z. In addition, it is noted tiiat since 
tiie vertical gradient is calculated using data from two different receivers witii two 
different ti-ansmitter positions, the noise in the two data sets will be incoherent. As a 
consequence, the noise level in the vertical gradient can be reduced. 

Although the horizontal field components of electric field and magnetic flux 
density are botii TM and TE dependant, the combinations shown in Equations 7 and 8 
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depend only on the TM mode. For this reason the combined response data defined by 
Equations 7 and 8 are referred to as TM mode decomposition data. In particular, 
because the TM mode decomposition shown in Equation 7 is maximum for inline 
orientations (i.e. cos(<|>) dependence), it is referred to as inline TM mode 
decomposition data. Because the TM mode decomposition shown in Equation 8 is 
maximum for broadside orientations (i.e. sin(^) dependence), it is referred to as 
broadside TM mode decomposition data. There are two significant benefits which are 
provided by combining the electric and magnetic field data in this way. Firstly, 
because the TM mode decomposition data does not include any dependence on the TE 
mode, the TM mode decomposition data is much less sensitive to the airwave 
component which prevents conventional analysis methods from working well in 
shallow water. Secondly, because Equation 8 provides combined response data having 
solely a TM mode component for broadside orientations, response data which are 
sensitive to the presence of a subterranean hydrocarbon reservoir may be obtained for 
broadside, as well as inline, transmitter-receiver orientations. This greatly increases 
the amount of usefiil data which can be obtained during a survey. Between them, 
Equations 7 and 8 can provide combined response data for all transmitter-receiver 
orientations. 

Figure 8 shows in schematic vertical section a model background subterranean 
strata configuration. The configuration comprises a section of seafloor 6 beneath a 
depth h of seawater 4. The seawater has a resistivity of 0.3 Qm. Beneath the seafloor 6 
is a uniform half-space sedimentary structure with a resistivity of 1 Qm, the low 
resistivity being primarily due to aqueous saturation of pore spaces. This background 
subterranean strata configuration extends uniformly downwards for an infinite extent. 
Also indicated in Figure 8 are an HED transmitter 22, and a receiver 25 of the kind 
shown in Figure 1. The distance between the HED transmitter and the receiver (i.e. the 
separation or range) is r km. The receiver is aligned inline (i.e. <|) = 0) with the HED 
transmitter. 
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Figure 9A is a graph schematically showing the logarithm of the modelled 
radial electric field component amplitude, Logio(E), seen at the receiver 25 in 
response to the HED transmitter 22 broadcast signal as a function of separation, r, 
between the transmitter and the receiver. This is the field component previously used 
as the basis for the analysis of CSEM survey data and is shown here for comparison 
purposes. Curves are calculated for a number of different water depths (H = 1500 m, 
1000 m, 500 m, 200 m and 100 m) as indicated on the figure. The HED transmitter is 
driven by an AC drive signal at a fi-equency of 0.25 Hz and the electric fields are 
calculated per unit transmitter electric dipole moment. Figure 9A demonstrates how 
the radial component of the electric field given by Equation 1 becomes increasingly 
dominated by the airwave component of the transmitted signal in shallower water. For 
example, at a separation of around 9000 m, the calculated radial electric field is 
approximately 300-times greater in a water depth of 100 m than in a water depth of 
1500 m. This is due to the increased relative contribution of the airwave component. 
Even at separations of only around 2000 m, the increased airwave contribution seen 
with a water depth of 100 m leads to radial electric fields which are around ten-times 
greater than those seen in deeper waters. 

Figure 9B is a graph schematically showing the phase, x, relative to the HED 
transmitter AC drive signal, of the modelled radial electric field components plotted in 
Figure 9A. It is apparent Grom Figure 9B that with a finite water depth there is little 
advance in phase with increasing separation once the airwave component begins to 
dominate, for example at beyond around r = 2000 m for h = 1 00 m. This is because a 
dominant component of the signal is travelling through the non-conducting air. 

Figure 1 OA is a graph schematically showing the logarithm of the modelled 
inline TM mode decomposition given by Equation 7 seen at the receiver 25 in 
response to the HED transmitter 22 broadcast signal multiplied by transmitter-receiver 
separation, r, as a fiinction of this separation. The multiplication by r provides for an 
equivalent electric field parameterisation of the TM mode decomposition. As with 
Figure 9A, curves are calculated for a number of different water depths h. The HED 
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transmitter is again driven by an AC drive signal at a frequency of 0.25 Hz and the 
TM mode decomposition calculated per unit transmitter electric dipole moment. It is 
clear from Figure lOA that, unlike Figure 9A, there is little difference between the 
curves for the different water depths. This reflects the fact that the TM mode 
decomposition does not include a TE mode dependence which is the mode which 
contributes most to the airwave component. 

Figure 1 OB is a graph schematically showing the phase, Xj relative to the HED 
transmitter AC drive signal, of the modelled TM mode decomposition plotted in 
Figure lOA. It is apparent from Figure lOB that the phase advances steadily with 
increasing separation for all water depths. This again demonstrates the insensitivity of 
the TM mode decomposition given by Equation 7 to the airwave component in 
shallow water. 

Figures IOC and lOD are similar to and will be understood from Figures lOA 
and lOB respectively. However, whereas Figures lOA and lOB show the behaviour of 
the modelled inline TM mode decomposition given by Equation 7 for the subterranean 
strata configuration shown in Figure 8, Figures IOC and lOD show the behaviour of 
the modelled broadside TM mode decomposition given by Equation 8. This is for the 
same subterranean strata configuration shown in Figure 8, but in which the receiver is 
now aligned broadside (i.e. ^ = 90) to the HED transmitter. 

This demonstration that both the inline and broadside TM mode 
decompositions are not significantly effected by the airwave component in shallow 
water is made using the model background subterranean strata configuration shown in 
Figure 8. This model does not contain a hydrocarbon reservoir. It is, therefore, 
important to show that the TM mode decompositions are sensitive to the presence of a 
hydrocarbon reservoir if they are to be of practical use. 

Figure 1 1 shows in schematic vertical section a model hydrocarbon-reservoir 
subterranean strata configuration. A section of seafloor 6 lies beneath a 1 00 m depth 
of seawater 4 which has a resistivity of 0.3 Qm. The strata configuration beneath the 
seafloor 6 comprises a 1000 m thick overburden layer 8, representing sediments. 



-29- 



arranged above a hydrocarbon reservoir 12. The overburden layer 8 has a resistivity of 
1 Qm, again, primarily due to aqueous saturation of pore spaces. The hydrocarbon 
reservoir 12 is 100 m thick, and has a resistivity of p Qm. This resistivity will 
typically be greater than that of the surrounding layers due to the presence of non- 
conducting hydrocarbon within pore spaces. Below the hydrocarbon reservoir 12 is a 
sedimentary underburden layer 9, which, as for the overburden layer, has a resistivity 
of 1 Qm. The underburden layer extends downwardly for an effectively infinite 
extent. Accordingly, except for the presence or absence of the hydrocarbon reservoir 
12, the hydrocarbon-reservoir subterranean strata configuration of Figure 11 is 
identical to the background subterranean strata configuration of Figure 8 for the case h 
= 100 m. An HED transmitter 22 and a receiver 25 are again shown. 

Figure 12A is a graph schematically showing the logarithm of the modelled 
radial electric field component amplitude, Logio(E), seen at the receiver 25 in 
response to the HED transmitter 22 broadcast signal as a function of separation, r, 
between the transmitter and the receiver with the hydrocarbon-reservoir subterranean 
reservoir shown in Figure 1 1 . This field component is again shown for comparison 
purposes. Curves are calculated for a number of different resistivities p for the 
hydrocarbon reservoir (p = 1 Qm (i.e. effectively no detectable reservoir), 10 Qm, 20 
Qm, 50 Qm and lOOQm) as indicated on the figure. The HED transmitter is again 
driven by an AC drive signal at a firequency of 0.25 Hz and the electric fields are 
calculated per unit transmitter electric dipole moment. The curves shown in Figure 
12A are all very similar to one another, even though there is a wide range of 
hydrocarbon resistivities. This is because with a water depth of only 1 00 m, the radial 
electric field component is dominated by the airwave component of the TE mode and 
cannot be used to properly identify the presence or not of a hydrocarbon reservoir. 

Figure 12B is a graph schematically showing the phase, x> relative to the HED 
transmitter AC drive signal, of the modelled radial electric field components plotted in 
Figure 12 A. It is apparent firom Figure 12B that there is little advance in phase with 
increasing separation for all hydrocarbon reservoir resistivities. This is again because 
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a dominant component of the transmitted signal is travelling through the non- 
conducting air. 

Figure 13A is a graph schematically showing the logarithm of the modelled 
inline TM mode decomposition seen at the receiver 25 in response to the HED 
5 transmitter 22 broadcast signal multiplied by transmitter-receiver separation as a 
function of this separation for the hydrocarbon-reservoir subterranean strata 
configuration of Figure 11 . As before, the multiplication by r provides an equivalent 
electric field parameterisation of the TM mode decomposition. As with Figure 12 A, 
curves are calculated for a number of different hydrocarbon reservoir resistivities. The 

10 HED transmitter is again driven by an AC drive signal at a firequency of 0.25 Hz and 
the TM mode decomposition calculated per unit transmitter electric dipole moment. It 
is clear from Figure 13A that, unlike the airwave dominated radial electric field curves 
of Figure 12 A, there is a strong dependence in the calculated response of the inline 
TM mode decomposition on the resistivity of the hydrocarbon reservoir, even though 

15 the seawater depth is only 100 m. For a hydrocarbon-reservoir resistivity of p = 100 
Qm, the inline TM mode decomposition signal is around 100-times greater at a 
separation of r = 11000 m than for the case p = 1 Qm (i.e. effectively no detectable 
hydrocarbon reservoir). This clearly demonstrates the sensitivity of the inline TM 
mode decomposition to the presence or not of a hydrocarbon reservoir. 

20 Figure 1 3B i s a graph schematically showing the phase, relative to the HED 

transmitter AC drive signal, of the modelled TM mode decomposition plotted in 
Figure 13 A. It is apparent fi-om Figure 13B that the phase advances at different rates 
for different hydrocarbon-reservoir resistivities. This again demonstrates the 
sensitivity of the inline TM mode decomposition given by Equation 7 to the presence 

25 of a hydrocarbon reservoir. 

Figures 13C and 13D are similar to and will be understood from Figures 13A 
and 13B respectively. However, whereas Figures 13 A and 13B show the behaviour of 
the modelled inline TM mode decomposition given by Equation 7 for the subterranean 
strata configuration shown in Figure 1 1, Figures 13C and 13D show the behaviour of 
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the modelled broadside TM mode decomposition given by Equation 8 multiplied by 
the transmitter-receiver separation in order to provide an equivalent electric field. This 
is for the same subterranean strata configuration as shown in Figure 1 1 , but in which 
the receiver is now aligned broadside (i.e. <|) = 90) to the HED transmitter. 

It can be seen fiom Figures 13C and 13D that there is also a strong dependence 
in the calculated response of the broadside TM mode decomposition on the resistivity 
of the hydrocarbon reservoir. For a hydrocarbon-reservoir resistivity of 100 Qm, the 
broadside TM mode decomposition signal is aroimd 300-times greater at a separation 
of r = 1 1000 m than for the case p = 1 ilm. This clearly demonstrates the sensitivity of 
the broadside TM mode decomposition to the presence or not of a hydrocarbon 
reservoir. 

Figure 14A is a graph schematically showing the ratio P of the curves plotted 
in Figures 12 A (radial electric field component), 13A (inline TM mode 
decomposition) and 13C (broadside TM mode decomposition) for a hydrocarbon 
reservoir of resistivity p = 100 Qm to the corresponding curves for which there is no 
detectable hydrocarbon reservoir (i.e. p — 1 ^2m). The curves are marked Er, TM 
and TM* " for the radial electric field, inline TM mode decomposition and broadside 
TM mode decomposition respectively. Figure 14A demonstrates the sensitivity of 
both the inline and broadside TM mode decompositions to the presence of the 
hydrocarbon reservoir as a fimction of separation r. This is apparent firom the large 
diversions firom unity for these curves. As noted above, at a separation of r = 1 1000 m, 
the inline TM mode decomposition is around 1 00-times greater, and the broadside TM 
mode decomposition is around 300-times greater with a p = 100 Qm hydrocarbon 
reservoir than when there is no detectable hydrocarbon reservoir (i.e. p = 1 Qm). The 
insensitivity of the radial electric field component to the presence of the hydrocarbon 
reservoir (due to the airwave component dominating the signal) is also clear. 

Figure 14B is a graph schematically showing the difference in phase Ax 
between the curves plotted in Figures 12B (radial electric field component), 13B 
(inline TM mode decomposition) and 13D (broadside TM mode decomposition) for a 
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hydrocarbon reservoir of resistivity p = 100 Qm and the respectively corresponding 
curves for which there is no detectable hydrocarbon reservoir (i.e. p = 1 Qm). The 
curves are marked Er, TM* ^ ^ and TM* ^ ^ for the radial electric field, inline TM 
mode decomposition and broadside TM mode decomposition respectively. Figure 14B 
again demonstrates the sensitivity of both the inline and broadside TM mode 
decompositions to the presence of the hydrocarbon reservoir as a function of 
separation r. This is apparent firom the progressive increase in the absolute value of 
Ax. The relative insensitivity of the radial electric field component to the presence of 
the hydrocarbon reservoir is again seen. 

Figure 14C is a graph schematically showing the same curves as shown in 
Figure 14A, but on a magnified scale. This more clearly demonstrates the sensitivity 
of the inline and broadside TM mode decompositions to the presence of the 
hydrocarbon reservoir at smaller separations r than can be seen on the scale of Figure 
14A. 

Curves of the kind shown in Figures 12 to 14 which are derived fi"om actual 
CSEM response data can be further analysed using standard techniques, for example, 
geophysical inversion, to produce subterranean resistivity maps of the area being 
surveyed. These analysis techniques can be broadly similar to techniques previously 
used in deep water surveys for electric field data of the kind shown in Figure 12A for 
conventional CSEM survey data analysis techniques for surveys, for example. 

Because in practice, subterranean strata configurations are generally not as 
simple as those used in the model surveys described above, it is sometimes difficult to 
identify directly from curves of the type shown in Figures 13 A, 13B, 13C and 13D 
obtained from real surveys whether the curves contain features indicative of a buried 
hydrocarbon reservoir or merely features relating to local larger scale background 
structures. In particular the kind of TM mode decomposition data seen with a thin 
resistive hydrocarbon reservoir embedded in a uniform resistivity background can be 
similar to that seen in a subterranean strata configuration comprising layers of 
increasing resistivity with depth. This kind of increasing-resistivity structure is a 
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feature of some submarine sedimentary basins, for example, and can arise due to the 
progressive expulsion of conductive pore fluids with increasing depths by a rising 
overburden pressure. Accordingly knowledge of the large scale background structure 
of the subterranean strata in the area from which survey data are being analysed is 
often helpful in order to determine reliably whether features in TM mode 
decomposition data are caused by a buried hydrocarbon layer or whether they are 
caused by large scale background structures. 

Figures 15A and 15B show two subterranean strata model configurations used 
to show the difficulty in distinguishing between a thin resistive hydrocarbon reservoir 
(Figure 15 A) and a steadily increasing resistivity with increasing depth (Figure 15B). 
Figure 1 5A shows a hydrocarbon-reservoir subterranean model configuration which is 
similar to that of Figure 1 1 for the case where the hydrocarbon reservoir resistivity p = 
100 Qm. However, the model subterranean strata configuration of Figure 15 A 
includes an infinite depth of seawater, as opposed to the 100 m depth of seawater of 
Figure 11. In the increasing-resistivity subterranean strata configuration model of 
Figure 15B, a section of seafloor 6 lies beneath an infinite depth of seawater 4. The 
strata beneath the seafloor 6 comprise a series of sedimentary layers of increasing 
resistivity. A first layer 10 has a uniform resistivity of 1 Qm and a thickness of 400 m. 
A second layer 13 has a uniform resistivity of 5 Qm and a thickness of 1000 m. 
Beneath the second layer 13 is a third layer 15 which has a resistivity of 10 Qm and 
extends downwardly for an infinite extent. An HED transmitter 22 and a receiver 25 
are also shown. 

Figure 16 is a graph showing modelled curves for the inline and broadside TM 
mode decomposition data which are similar to and will be understood from the TM 
mode decomposition curves shown in Figure 14 A, but which are calculated for the 
hydrocarbon-reservofr subterranean strata configuration shown in Figure 15A (solid 
lines) and for the increasing-resistivity subterranean strata configuration shown in 
Figure 15B (dotted lines). It is clear that the inline TM mode decomposition data 
calculated for the hydrocarbon reservoir model (solid line marked TM* ^ ^) are similar 
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to the inline TM mode decomposition data calculated for the increasing resistivity 
model (dotted line marked TM* " Likewise, the broadside TM mode decomposition 
data calculated for the hydrocarbon reservoir model (solid line marked TM* " are 
similar to the broadside TM mode decomposition data calculated for the increasing 
resistivity model (dotted line marked TM* " ^^). This demonstrates the ambiguity that 
can arise with TM mode decomposition data when attempting to distinguish between 
a subterranean strata configuration having a hydrocarbon reservoir and some other 
large scale subterranean strata configurations. 

Because of this possible ambiguity, analysis of survey data aimed at 
establishing whether a subterranean strata configuration contains a thin resistive 
hydrocarbon reservoir will normally involve combining electric and magnetic field 
data to generate TM mode decomposition data such as that defined by Equations 7 or 
8. These response data are sensitive to the presence of subterranean hydrocarbon 
reservoirs, even in shallow seawater. However, in addition, to determine reliably 
whether features of the TM mode decomposition data are indicative of a hydrocarbon 
reservoir or of the local background structure, it is necessary to determine how the TM 
mode decomposition data for a given subterranean strata configuration would appear 
if there were no hydrocarbon reservoir. 

This analysis step, generally referred to as normalisation, is usually done with 
the aid of background data. Background data are specific to the area being surveyed 
and can be obtained in a variety of ways. One way is to model the EM survey 
performed to obtain the TM mode decomposition data with a model background 
subterranean strata configuration. The background model should be as close a match 
as possible to the actual background structure in the area being surveyed. A 
comparison of the TM mode decomposition data with the background data provides 
difference data sensitive to the likely presence, extent and location of a subterranean 
hydrocarbon reservoir embedded within the background subterranean strata 
configuration. For example, if the TM mode decomposition data closely match the 
background data, there is unlikely to be a buried hydrocarbon layer. If, on the other 
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hand, there are differences, i.e. anomalies, in the TM mode decomposition data 
compared to the background data, for example, an increased receiver signal 
amplitude, this could be quantitatively assessed in terms of being indicative of a 
buried hydrocarbon reservoir. The variation in anomalies at different horizontal 
separations provides information on the depth and extent of a hydrocarbon reservoir 
For example, if differences between the TM mode decomposition data and the 
background data are only apparent at large transmitter-receiver horizontal separations, 
this is likely to indicate that the hydrocarbon reservoir is relatively deeply buried. 
Similarly, a discontinuity in TM mode decomposition data as a function of horizontal 
separation is likely to indicate a boundary or edge of a hydrocarbon reservoir at the 
location of the discontinuity. 

Suitable background models to use in generating background data can be 
obtained in several ways. 

One way of obtaining the information required to construct a suitable 
background model is with conventional MT electromagnetic surveying techniques. As 
noted above, these techniques are capable of providing information on large scale 
background resistivity structures, even though they unable to detect hydrocarbon 
reservoirs directly. 

Another way of obtaining the information required to construct a suitable 
background model is from CSEM survey data. As described above, it is the TE mode 
component of a transmitted signal which can provide information on background 
structure is a CSEM survey. TE mode response may be obtained from a combination 
of the same electric and magnetic field data combined to provide the TM mode 
decomposition data. 

Equation 9, shown in Figure 1 7A, defines a linear combination of the radial 
electric field data Er and the vertical gradient in the azimuthal magnetic field data B^. 
Equation 9 defines combined response data which include only a TE mode 
dependence and no TM mode dependence. 
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Equation 1 0, shown in Figure 1 7B, defines a similar linear combination of the 
azimuthal electric field data and the vertical gradient in the radial magnetic field 
data Bf. This combined response data again includes only a TE mode dependence. 

The vertical gradients in the magnetic field data may be obtained similarly to 
5 the vertical gradients in the electric field described above. 

The combined response data given by Equation 9 is referred to as inline TE 
mode decomposition data and the combined response data shown in Equation 10 is 
referred to as broadside TE mode decomposition data. 

Figure 1 8A is a graph which is similar to and will be understood fi"om Figure 
10 16. However, whereas Figure 16 plots data for the inline and broadside TM mode 
decompositions calculated for the model subterranean strata shown in Figures ISA 
and 1 5B, Figure 1 8 A plots data for the inline and broadside TE mode decompositions 
calculated for the same model subterranean strata configurations. Curves calculated 
for the hydrocarbon-reservoir subterranean strata configuration shown in Figure 1 5 A 
15 are shown as solid lines and curves for the increasing-resistivity subterranean strata 
configuration shown in Figure 1 5B are shown as dotted lines. It is clear that the inline 
TE mode decomposition data calculated for the hydrocarbon reservoir model (solid 
line marked TE* " ^) is very different to the inline TE mode decomposition data 
calculated for the increasing resistivity model (dotted line marked TE* " % Likewise, 
20 the broadside TE mode decomposition data calculated for the hydrocarbon reservoir 
model (solid line marked TE* ^ is very different to the broadside TE mode 
decomposition data calculated for the increasing resistivity model (dotted line marked 
TE* " This is because, as described above, the TE mode decomposition is 
insensitive to the hydrocarbon reservoir for both inline and broadside alignments. 
25 Figure 18B is a graph schematically showing the same curves as shown in 

Figure 1 8 A, but on a magnified scale. 

With conventional survey analysis methods relying on the amplitude of 
detected electric fields, the TE mode response is obtained by aligning the transmitter 



-37- 

and the receiver in a broadside orientation. It is not possible with conventional survey 
analysis methods to determine the TE mode response at inline orientations. 

Whilst in shallow water the use of TE mode decompositions in determining a 
background model is prone to the same difficulties associated with the airwave 
5 component described above, the impact of the airwave component can be reduced by 
employing relatively low fi^equency EM signals. Low frequency signals suffer less 
attenuation as they pass through the subterranean strata and so the airwave component 
is not so dominant in the EM fields induced at a receiver by an HED transmitter 
driven by a low frequency AC current. Because of this, low frequency signals are 

10 capable of providing information on large scale backgroimd resistivity structures 
needed to generate a background model. 

hi other cases, an area to be surveyed will already be very well characterized 
by previous surveying. For example, in a producing oilfield or oil province there is 
likely to be a wealth of existing seismic and well-log data. In these cases, background 

1 S models can be calculated from a rock formation model. The rock formation model can 
be created from the seismic data and then resistivities assigned to the various 
components in the rock structure using the resistivities obtained from well-log data. (If 
directly applicable well-log data are not available, it may be possible to estimate 
resistivity values by comparison with resistivity data from nearby wells in similar 

20 geological structures.) This technique for obtaining the information required to 
constmct a suitable background model will be especially suited to applications in 
existing oilfields, such as monitoring long term depletion of reserves. 

When monitoring depletion, it may be sufficient to directly compare TM mode 
decomposition data taken at different times, e.g. several weeks or months apart, 

25 without use of a rock formation model. In other words, the background data used is 
data from a previous similar survey. Differences in TM mode decomposition data 
taken at different times are indicative of changes in the hydrocarbon reservoir which 
have occurred between the times at which the data were taken. Because of this, this 
kind of comparison provides a usefiil monitoring tool. The TM mode decomposition 
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data taken at the earlier time thus effectively acts as background data for comparing 
with the TM mode decomposition data taken at the later time. 

In addition to the insensitivity of the TM mode decompositions to the airwave 
component, another major benefit is that TM mode decompositions may be made with 
5 both inline (Equation 7) and broadside (Equation 8) orientations. Furthermore, 
because the TM mode decompositions are not mixed mode at any orientation, the 
response of the subterranean strata to the TM mode can be properly determined for all 
transmitter-receiver orientations, the only effect being an effective reduction in the 
transmitter dipole amplitude. At transmitter-receiver orientations which are closer to 
10 inline than broadside, the TM mode decomposition according to Equation 7 provides 
the greatest signal. At transmitter-receiver orientations which are closer to broadside 
than inline, the TM mode decomposition according to Equation 8 provides the 
greatest signal. 

Similar considerations allow the response of the subterranean strata to the TE 
15 mode decomposition to be obtained for all transmitter-receiver orientations. This is 

useful where TE mode decomposition data are used for normalisation, for example 

when analysing survey data firom a deep-water survey. 

Because both TM and TE mode decomposition data can be obtained for all 

transmitter-receiver orientations, and not just the particular configuration required 
20 with conventional analysis methods, survey response data analysed according to the 

above described techniques can be obtained using much simpler tow paths than has 

previously been possible. 

Figure 19 is a schematic plan view showing an example layout of sixteen 

receivers 25 distributed across a section of seafloor 6 above a subterranean reservoir 
25 56 during a survey for which data are to be analysed according to the above described 

TM mode decomposition techniques. The reservoir 56 has a linear extent on the order 

of several km and its boundary is indicated by a heavy line 58. The reservoir is similar 

to the reservoir shown in Figure 2. The orientation of the subterranean reservoir is 

indicated by the cardinal compass points (marked N, E, S and W for North, East, 
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South and West respectively) indicated in the upper right of the figure. In this 
example, the receivers 25 are uniformly distributed in a square-grid pattern so as to 
approximately cover the reservoir 56. In performing a survey, an HED transmitter (not 
shown) starts from the location marked 'A' in Figure 19 and is towed, whilst 
5 broadcasting continuously as described above, along a path indicated by the broken 
line 60. The survey is completed when the transmitter reaches the location marked 
*B'. Electric and magnetic field data are continuously collected by the receivers 25 
throughout the towing process and the position of the HED transmitter relative to the 
network of receivers is also logged. 

10 During the towing process, each of the receivers 25 presents several different 

transmitter-receiver horizontal separations and orientations relative to the HED 
transmitter. Accordingly, by following the tow path marked, electric and magnetic 
field data are collected for many different transmitter-receiver horizontal separations 
along many different directions. These data can be combined to generate combined 

1 5 response data of the kind defined in Equation 7 or 8. The combined response data may 
be inverted to provide a thorough map of the subterranean strata configuration using 
appropriate geophysical inversion techniques. In this way a simple continuous towing 
of the transmitter can provide a detailed survey which covers the extent of the 
subterranean reservoir 56. In cases where vertical gradients in the fields are to be 

20 determined by towing a conventional HED transmitter at two different heights, the 
path shown in Figure 19 will be followed twice. Firstly with the HED transmitter at a 
first height and then subsequently with the transmitter at a second height. Although in 
general it will be preferable for the HED transmitter to be towed at a fairly constant 
height during each tow, conventional radar and/or echo location techniques, for 

25 example, can be used to record the height of the HED transmitter at each location 
during each tow. This means vertical gradients in the fields can be determined based 
on differences in the detected signals and differences in the measured heights at each 
location. 
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As described above, Figure 2 shows for comparison an example tow path 
which might be followed in performing a survey to obtain data to be analysed using 
previously known techniques. The tow path shown in Figure 2 is approximately twice 
as long as that shown in Figure 19. Furthermore, when performing an EM survey 
according Figure 1 9, useful response data can be collected when the HED transmitter 
is at all points along a tow path. This is because a suitable TM mode decomposition 
can be deteraiined for all transmitter-receiver orientations. 

The TM mode decomposition analysis technique thus allows a greater amount 
of useful data to be collected from a shorter tow path than the previous proposed end- 
on/broadside analysis method. The same is tme of the TE mode decompositions 
where these are used for background normahsation, for example in deep water surveys 
where the airwave effect on the TE mode decomposition is not important. In these 
cases TM and TE mode decomposition data can be obtained simultaneously and for 
the same transmitter-receiver separations and orientations. This has the further benefit 
of reducing any ambiguities arising where the subterranean strata is not truly one 
dimensional. 

Although the above example is based on a square receiver grid, it will be 
understood that a wide variety of receiver placements may be used. For example other 
high symmetry regular grids, such as triangular or rectangular, may be used. In 
addition irregular grids may be used that have no high level of symmetry. 

During surveying, transmitter and receiver positions can be determined using 
standard long and/or short and/or ultra-short baseline acoustic navigation systems and 
processed to determine the horizontal separation between the transmitter and receiver. 

It can be seen from Figure 19 that with the presently proposed analysis 
techniques, subterranean strata can be comprehensively surveyed with relatively few 
receivers and EM transmitter tows. This reduces the time taken, and hence cost, of an 
exploratory EM survey. 

It will be understood that whilst the above description describes a towed HED 
transmitter, the method would also be applicable in a fixed installation. For example. 
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the method could be used to monitor changes to a hydrocarbon reservoir from which 
hydrocarbon is being drawn. In such cases it will be appropriate to employ one (or 
more) HED transmitter(s) in fixed positions relative to a receiver array rather than 
undertake frequent towed surveys. The HED transmitter(s) could be anchored to the 
seafloor or suspended from an oil-rig platform, for example. In other examples, the 
HED transmitter(s) could be placed in a horizontal well or borehole, e.g. a 
geotechnical borehole. In the case of a producing oil field, the subterranean structures 
are Ukely to be well known aheady from prior geophysical surveys and drilling 
results. Indeed, prior geophysical and geological information from the oil field can be 
used to construct a backgroimd model as described above. 

Finally, it will be understood that the invention is equally applicable to 
surveying of freshwater, for example large lakes or estuaries, so that references to 
seafloor, seawater etc. should not be regarded as limiting and should be interpreted as 
covering lakebed, riverbed etc. Indeed the applicability of the invention to shallow 
water makes it ideal for surveying shallow lakes. 
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CLAIMS 

1 . A method of analysing results from an electromagnetic survey of an area that 
is thought or known to contain a subterranean hydrocarbon reservoir, comprising: 

5 providing electric field data and magnetic field data obtained by at least one 

receiver from at least one horizontal electric dipole (HED) transmitter; 
determining a vertical gradient in the electric field data; and 
combining the vertical gradient in the electric field data with the magnetic 
field data to generate combined response data. 

10 

2. A method of analysing results from an electromagnetic survey according to 
claim 1 , wherein the electric field data include a horizontal component of electric field 
resolved along a first direction and the magnetic field data include a horizontal 
component of magnetic field data resolved along a second direction, the first and 

1 5 second directions being different 

3. A method of analysing resuhs from an electromagnetic survey according to 
claim 2, wherein the first and second directions are orthogonal to one another. 

20 4. A method of analysing results fix)m an electromagnetic survey according to 
claim 2 or claim 3, wherein the first direction is parallel to a line connecting the HED 
transmitter to the receiver. 

5. A method of analysing results from an electromagnetic survey according to 
25 claim 2 or claim 3, wherein the first direction is perpendicular to a line connecting the 
HED transmitter to the receiver. 
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6. A method of analysing results from an electromagnetic survey according to 
any of claims 1 to 5, wherein the vertical gradient in the electric field data is 
determined by comparing electric field data detected at different heights. 

5 7. A method of analysing results fix^m an electromagnetic survey according to 
any of claims 1 to 5, wherein the vertical gradient in the electric field data is 
determined by comparing the electric field data and data simulated using a 
background model. 

10 8. A method of analysing results from an electromagnetic survey according to 
claim 7, wherein the data simulated using a background model provide a boundary 
condition for the electric field data. 

9. A method of analysing results from an electromagnetic survey according to 
15 any of claims 1 to 5, wherein the vertical gradient in the electric field data at a first 

receiver is determined by comparing electric field data from the first receiver when the 
transmitter is above a second receiver with electric field data from the second receiver 
when the transmitter is above the first receiver, and applying a predetermined 
adjustment to the electric field data from second receiver. 

20 

10. A method of analysing results from an electromagnetic survey according to 
any of claims 1 to 5, wherein the vertical gradient in the electric field data is 
determined by comparing electric field data detected from a transmitter at different 
heights. 

25 

11. A method of analysing results from an electromagnetic survey according to 
any of claims 1 to 1 0, further comprising: 

providing background data specific to the area being surveyed; and 
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comparing the combined response data with the background data to obtain 
difference data sensitive to the presence of a subterranean hydrocarbon reservoir. 

12. A method of analysing results from an electromagnetic survey according to 
5 claim 1 1, wherein the background data are obtained by determining a vertical gradient 

in the magnetic field data and combining the vertical gradient in the magnetic field 
data with the electric field data. 

13. A method of analysing results from an electromagnetic survey according to 
10 claim 12, wherein the vertical gradient in the magnetic field data is determined by 

comparing magnetic field data detected at different heights. 

14. A method of analysing results from an electromagnetic survey according to 
claim 12, wherein the vertical gradient in the magnetic field data is determined by 

IS comparing the magnetic field data and data simulated using a background model. 

15. A method of analysing results from an electromagnetic survey according to 
claim 1 4, wherein the data simulated using a background model provide a boundary 
condition for the magnetic field data. 

20 

16. A method of analysing results from an electromagnetic survey according to 
claims 12, wherein the vertical gradient in the magnetic field data at a first receiver is 
determined by comparing magnetic field data from the first receiver when the 
transmitter is above a second receiver with magnetic field data from the second 

25 receiver when the transmitter is above the first receiver, and applying a predetemiined 
adjustment to the magnetic field data from second receiver. 
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17. A method of analysing results fiom an electromagnetic survey according to 
claims 12, wherein the vertical gradient in the magnetic field data is determined by 
comparing magnetic field data detected fi*om a transmitter at different heights. 

5 18. A method of analysing results firom an electromagnetic survey according to 
claim 11, wherein the background data are obtained firom a controlled source 
electromagnetic survey. 

19. A method of analysing results from an electromagnetic survey according to 
10 claim 11, wherein the background data are obtained from a magneto-telluric 

electromagnetic survey. 

20. A method of analysing results from an electromagnetic survey according to 
claim 1 1 , wherein the background data are fiirther combined response data obtained 

1 5 from another electromagnetic survey of the area performed at a different time. 

21. A method of analysing results from an electromagnetic survey according to 
claim 1 1 , wherein the background data are calculated from a rock formation model. 

20 22. A method of analysing results firom an electromagnetic survey according to 
claim 21, wherein the rock formation model is derived from a combination of 
geological data and resistivity data, 

23. A method of analysing results from an electromagnetic survey according to 
25 claim 22, wherein the geological data are from seismological surveying. 

24. A method of analysing results from an electromagnetic survey according to 
claim 22 or 23, wherein the resistivity data are from well logging. 
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25. A method of analysing results from an electromagnetic survey according to 
any of claims 1 to 24, wherein difference data are obtained as a fimction of position 
within the area. 

5 26. A computer program product bearing machine readable instructions for 
implementing a method of analysing results from an electromagnetic survey according 
to any of claims 1 to 25. 

27. A computer apparatus loaded with machine readable instructions for 
10 implementing the method of analysing results from an electromagnetic survey 

according to any of claims 1 to 25. 

28. A method of planning an electromagnetic survey of an area that is thought or 
known to contain a subterranean hydrocarbon reservoir, comprising: 

15 creating a model of the area to be surveyed including a rock formation 

containing a postulated hydrocarbon reservoir, and a body of water above the rock 
formation; 

setting values for water depth, depth of the postulated hydrocarbon reservoir, 
and resistivity structure of the rock formation; and 
20 performing a simulation of an electromagnetic survey in the model of the 

survey area by calculating electric field data and magnetic field data obtained by at 
least one simulated receiver detecting signals from at least one simulated horizontal 
electric dipole (HED) transmitter; 

determining a vertical gradient in the electric field data; and 
25 combining the vertical gradient in the electric field data with the magnetic 

field data to generate combined response data. 

29. A method of planning an electromagnetic survey according to claim 28, further 
comprising: 
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adjusting the model to remove the postulated hydrocarbon reservoir; and 
repeating the simulation to obtain background data for comparison with the 
combined response data. 

30. A method of planning an electromagnetic survey according to claim 28 or 
claim 29, the method further comprising: 

repeating the simulation for a number of transmitter-receiver horizontal 
separations and frequencies of transmitter signal in order to select optimum surveying 
conditions in terms of transmitter-receiver horizontal separations and frequencies for 
probing the hydrocarbon reservoir. 

31. A computer program product bearing machine readable instructions for 
implementing the method of planning an electromagnetic survey according to any of 
claims 28 to 30. 

32. A computer apparatus loaded with machine readable instructions for 
implementing the method of planning an electromagnetic survey according to any of 
claims 28 to 30. 

33. An electromagnetic survey method applied to a survey area that is thought or 
known to contain a subterranean hydrocarbon reservoir, comprising: 

providing at least one transmitter and at least one detector for transmission and 
detection of electromagnetic signals; and 

obtaining data with transmission and/or detection at a plurality of different 
heights over the survey area, so that the data allow comparison of electromagnetic 
signals transmitted and/or received at different vertical displacements. 
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34. An electromagnetic sxu^ey method according to claim 33, wherein the 
transmission and/or detection at a plurality of different heights comprises detection at 
a plurality of different heights. 

5 35. An electromagnetic survey method according to claim 34, wherein the 
detection at a plurality of different heights is made simultaneously by a corresponding 
plurality of detectors. 

36. An electromagnetic survey method according to claim 33, wherein the 
10 transmission and/or detection at a plurality of different heights comprises transmission 

at a plurality of different heights. 

37. An electromagnetic survey method according to claim 36, wherein the 
transmission at a plurality of different heights is made simultaneously by a 

1 5 corresponding pliurality of transmitters. 

38. An electromagnetic survey method according to claim 36, wherein the 
transmission at a plurality of different heights is made at different times by a single 
transmitter. 

20 

39. An electromagnetic survey method according to any of claims 33 to 38, 
wherein data are obtained as function of position over the survey area. 

40. An electromagnetic receiver for use in an electromagnetic survey of an area 
25 that is thought or known to contain a subterranean hydrocarbon reservoir, the receiver, 

when normally deployed, being operable to measure electric fields at two or more 
different heights such that a vertical gradient in electric field may subsequently be 
determined. 
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41. An electromagnetic receiver according to claim 40, wherein the receiver is 
further operable to measure magnetic fields at two or more different heights such that 
a vertical gradient in magnetic field may subsequently be determined. 

5 42. An electromagnetic receiver according to claim 40 or 41, wherein the receiver, 
when normally deployed, comprises first and second pairs of horizontally aligned 
electric dipole detectors extending in different horizontal directions, the pairs being 
positioned at different heights. 

10 43. An electromagnetic receiver according to claim 41 or 42, wherein the receiver, 
when normally deployed, comprises first and second pairs of magnetic field detectors 
for detecting magnetic fields in different horizontal directions, the pairs being 
positioned at different heights. 

15 44. An electromagnetic source for use in an electromagnetic survey of an area that 
is thought or known to contain a subterranean hydrocarbon reservoir, the source 
comprising first and second transmitters which in normal use are disposed at different 
heights. 

20 45. A method of analysing results from an electromagnetic survey substantially as 
hereinbefore described with reference to Figures 4 to 19 of the accompanying 
drawings. 

46. A method of planning an electromagnetic survey substantially as hereinbefore 
25 described with reference to Figures 4 to 19 of the accompanying drawings. 

47. A method of electromagnetic surveying substantially as hereinbefore described 
with reference to Figures 6A and 6B of the accompanying drawings. 
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48, An electromagnetic receiver substantially as hereinbefore described with 
reference to Figure 6A of the accompanying drawings. 

49. An electromagnetic transmitter substantially as hereinbefore described with 
5 reference to Figure 6B of the accompanying drawings. 



Amendments to the claims have been filed as follows 



CLAIMS 

1 . A method of analysing results from an electromagnetic survey of an area that 
is thought or known to contain a subterranean hydrocarbon reservoir, comprising: 

providing electric field data and magnetic field data obtained by at least one 
receiver from at least one horizontal electric dipole (HED) transmitter; 
determining a vertical gradient in the electric field data; and 
combining the vertical gradient in the electric field data with the magnetic 
field data to generate combined response data. 

2. A method of analysing results from an electromagnetic survey according to 
claim 1 , wherein the electric field data include a horizontal component of electric field 
resolved along a first direction and the magnetic field data include a horizontal 
component of magnetic field data resolved along a second direction, the first and 
second directions being different. 

3. A method of analysing results firom an electromagnetic survey according to 
claim 2, wherein the first and second directions are orthogonal to one another. 

4. A method of analysing results fi-om an electromagnetic survey according to 
claim 2 or claim 3, wherein the first direction is parallel to a line connecting the HED 
transmitter to the receiver. 

5. A method of analysing results fi-om an electromagnetic survey according to 
claim 2 or claim 3, wherein the first direction is perpendicular to a line connecting the 
HED transmitter to the receiver. 



6. A method of analysing results from an electromagnetic survey according to 
any of claims 1 to 5, wherein the vertical gradient in the electric field data is 
determined by comparing electric field data detected at different heights. 

7. A method of analysing results from an electromagnetic survey according to 
any of claims 1 to 5, wherein the vertical gradient in the electric field data is 
deteraiined by comparing the electric field data and data simulated using a 
background model. 

8. A method of analysing results from an electromagnetic survey according to 
claim 7, wherein the data simulated using a background model provide a boundary 
condition for the electric field data. 

9. A method of analysing results from an electromagnetic survey according to 
any of claims 1 to 5, wherein the vertical gradient in the electric field data at a first 
receiver is determined by comparing electric field data from the first receiver when the 
transmitter is above a second receiver with electric field data from the second receiver 
when the transmitter is above the first receiver, and applying a predetermined 
adjustment to the electric field data fi"om second receiver. 

1 0. A method of analysing results from an electromagnetic survey according to 
any of claims I to 5, wherein the vertical gradient in the electric field data is 
determined by comparing electric field data detected fi-om a transmitter at different 
heights. 

11. A method of analysing results from an electromagnetic survey according to 
any of claims 1 to 10, further comprising: 

providing background data specific to the area being surveyed; and 



comparing the combined response data with the background data to obtain 
difference data sensitive to the presence of a subterranean hydrocarbon reservoir. 



12. A method of analysing results from an electromagnetic survey according to 
claim 1 1, wherein the background data are obtained by determining a vertical gradient 
m the magnetic field data and conibinii^g the vertical gradient in the magnetic field 
data with the electric field data. 



13. A method of analysing results from an electromagnetic survey according to 
claim 12, wherein the vertical gradient in the magnetic field data is determined by 
comparing magnetic field data detected at different heights. 

14. A method of analysing results from an electromagnetic survey according to 
claim 12, wherein the vertical gradient in the magnetic field data is determined by 
comparing the magnetic field data and data simulated using a background model. 

15. A method of analysing results from an electromagnetic survey according to 
claim 14, wherein the data simulated using a background model provide a boundary 
condition for the magnetic field data. 

16. A method of analysing results from an electromagnetic survey according to 
claims 12, wherein the vertical gradient in the magnetic field data at a first receiver is 
determined by comparing magnetic field data from the first receiver when the 
transmitter is above a second receiver with magnetic field data from the second 
receiver when the transmitter is above the first receiver, and applying a predetermined 
adjustment to the magnetic field data from second receiver. 



17. A method of analysing results from an electromagnetic survey according to 
claims 12, wherein the vertical gradient in the magnetic field data is determined by 
comparing magnetic field data detected from a transmitter at different heights. 

5 18. A method of analysing results from an electromagnetic survey according to 
claim 11, wherein the background data are obtained from a controlled source 
electromagnetic survey. 

19. A method of analysing results from an electromagnetic survey according to 
10 claim 11, wherein the background data are obtained from a magneto-telluric 

electromagnetic survey, 

20. A method of analysing results from an electromagnetic survey according to 
claim 1 1 , wherein the background data are further combined response data obtained 

15 from another electromagnetic survey of the area performed at a different time. 

21. A method of analysing results from an electromagnetic survey according to 
claim 1 1 , wherein the background data are calculated from a rock formation model. 

20 22. A method of analysing results from an electromagnetic survey according to 
claim 21, wherein the rock formation model is derived from a combination of 
geological data and resistivity data. 

23. A method of analysing results from an electromagnetic survey according to 
25 claim 22, wherein the geological data are from seismological surveying. 

24. A method of analysing results from an electromagnetic survey according to 
claim 22 or 23, wherein the resistivity data are from well logging. 



25. A method of analysing results from an electromagnetic survey according to 
any of claims 1 to 24, wherein difference data are obtained as a function of position 
within the area. 

26. A computer program product bearing machine readable instructions for 
implementing a method of analysing results from an electromagnetic sun^ey according 
to any of claims 1 to 25. 

27. A computer apparatus loaded with machine readable instructions for 
implementing the method of analysing results from an electromagnetic survey 
according to any of claims 1 to 25. 

28. A method of planning an electromagnetic survey of an area that is thought or 
known to contain a subterranean hydrocarbon reservoir, comprising: 

creating a model of the area to be surveyed including a rock formation 
containing a postulated hydrocarbon reservoir, and a body of water above the rock 
formation; 

setting values for water depth, depth of the postulated hydrocarbon reservoir, 
and resistivity structure of the rock formation; and 

performing a simulation of an electromagnetic survey in the model of the 
survey area by calculating electric field data and magnetic field data obtained by at 
least one simulated receiver detecting signals from at least one simulated horizontal 
electric dipole (HED) transmitter; 

determining a vertical gradient in the electric field data; and 

combining the vertical gradient in the electric field data with the magnetic 
field data to generate combined response data. 



29. A method of planning an electromagnetic survey according to claim 28, fiirther 
comprising: 



adjusting the model to remove the postulated hydrocarbon reservoir; and 
repeating the simulation to obtain background data for comparison with the 
combined response data. 

30. A method of planning an electromagnetic survey according to claim 28 or 
claim 29, the method further comprising: 

repeating the simulation for a number of transmitter-receiver horizontal 
separations and frequencies of transmitter signal in order to select optimum surveying 
conditions in terms of transmitter-receiver horizontal separations and frequencies for 
probing the hydrocarbon reservoir. 

31. A computer program product bearing machine readable instructions for 
implementing the method of planning an electromagnetic survey according to any of 
claims 28 to 30. 

32. A computer apparatus loaded with machine readable instructions for 
implementing the method of planning an electromagnetic survey according to any of 
claims 28 to 30. 

33. An electromagnetic survey method applied to a survey area that is thought or 
known to contain a subterranean hydrocarbon reservoir, the survey area comprising 
subterranean strata beneath a seafloor, the method comprising: 

providing at least one transmitter and at least one detector for transmission and 
detection of electromagnetic signals; and 

obtaining data with transmission and/or detection at a plurality of different 
heights above the seafloor over the survey area, so that the data allow comparison of 
electromagnetic signals transmitted and/or received at different vertical displacements. 



34. An electromagnetic survey method according to claim 33, wherein the 
transmission and/or detection at a plurality of different heights comprises detection at 
a plurality of different heights. 

35. An electromagnetic survey method according to claim 34, wherein the 
detection at a plurality of different heights is made simultaneously by a corresponding 
plurality of detectors. 

36. An electromagnetic survey method according to claim 33, wherein the 
transmission and/or detection at a plurality of different heights comprises transmission 
at a plurality of different heights. 

37. An electromagnetic survey method according to claim 36, wherein the 
transmission at a plurality of different heights is made simultaneously by a 
corresponding plurality of transmitters. 

38. An electromagnetic survey method according to claim 36, wherein the 
transmission at a plurality of different heights is made at different times by a single 
transmitter. 

39. An electromagnetic survey method according to any of claims 33 to 38, 
wherein data are obtained as function of position over the survey area. 

40. An electromagnetic receiver for use in an electromagnetic survey of an area 
that is thought or known to contain a subterranean hydrocarbon reservoir, the area 
comprising subterranean strata beneath a seafloor, the receiver, when normally 
deployed, being operable to measure electric fields at two or more different heights 
above the seafloor such that a vertical gradient in electric field may subsequently be 
determined. 



41. An electromagnetic receiver according to claim 40, wherein the receiver is 
further operable to measure magnetic fields at two or more different heights above the 
seafloor such that a vertical gradient in magnetic field may subsequently be 
determined. 

42. An electromagnetic receiver according to claim 40 or 41, wherein the receiver, 
when normally deployed, comprises first and second pairs of horizontally aligned 
electric dipole detectors extending in different horizontal directions, the pairs being 
positioned at different heights. 

43. An electromagnetic receiver according to claim 41 or 42, wherein the receiver, 
when normally deployed, comprises first and second pairs of magnetic field detectors 
for detecting magnetic fields in different horizontal directions, the pairs being 
positioned at different heights. 

44. An electromagnetic source for use in an electromagnetic survey of an area that 
is thought or known to contain a subterranean hydrocarbon reservoir, the area 
comprising subterranean strata beneath a seafloor, the source comprising first and 
second transmitters which in normal use are disposed at different heights above the 
seafloor. 

45. A method of analysing results from an electromagnetic survey substantially as 
hereinbefore described with reference to Figures 4 to 19 of the accompanying 
drawings. 

46. A method of planning an electromagnetic survey substantially as hereinbefore 
described with reference to Figures 4 to 1 9 of the accompanying drawings. 

47. A method of electromagnetic surveying substantially as hereinbefore described 
with reference to Figures 6A and 6B of the accompanying drawings. 



48. An electromagnetic receiver substantially as hereinbefore described 
reference to Figure 6A of the accompanying drawings. 

49. An electromagnetic transmitter substantially as hereinbefore described 
5 reference to Figure 6B of the accompanying drawmgs. 
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